Current Flow in Random Resistor Networks: 
The Role of Percolation in Weak and Strong Disorder 
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We study the current flow paths between two edges in a random resistor network on a L x L square 
lattice. Each resistor has resistance e ax , where x is a uniformly-distributed random variable and a 
controls the broadness of the distribution. We find (a) the scaled variable u = L/a v , where v is the 
percolation connectedness exponent, fully determines the distribution of the current path length t 
for all values of u. For u 2> 1, the behavior corresponds to the weak disorder limit and t scales as 
I ~ L, while for u <C 1, the behavior corresponds to the strong disorder limit with I ~ L d °& , where 
rfopt = 1-22 ± 0.01 is the optimal path exponent, (b) In the weak disorder regime, there is a length 
scale £ ~ a", below which strong disorder and critical percolation characterize the current path. 
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Transport in disordered media is a classic problem in 
statistical physics which attracts much attention due to 
its broad range of applications. Examples include flow 
through porous material, oil production, and conductiv- 
ity of semiconducting materials or metal-insulator mix- 
tures [IlSHSEIiHIlllIIS These Problems have 
been studied using a random resistor network model with 
bonds that have a resistance chosen from a probability 
distribution mimicking the nature of the physical prob- 
lem under consideration. Among the different classes of 
disorder distributions used, the most common is percola- 
tion disorder, in which the resistance of a bond is either 
1 or oo [Tlj . Gaussian distributions and power law dis- 
tributions have also been studied extensively |l2|, LL3| ■ 

Here, we stud y a random resistor network with expo- 
nential disorder [1J| . We consider the two opposite edges 
of a L x L square lattice as source A and sink B. Each 
bond connecting adjacent nodes i and j corresponds to 
one resistor, whose resistance r%j is given by jlll2ll.'ll4Lll. r j| 

rij =e ax «, (1) 

where a controls the disorder strength and Xij is a ran- 
dom number taken from a uniform distribution € 
[0,1]. Recent experiments show that for quenched con- 
densed granular Ni thin films, the conductivity is well 
described by exponential disorder with large a Q. Ex- 
ponential disorder enables us to understand the magne- 
toresistance phenomenon that out of 10 9 grains, only a 
few govern the electric conductivity §|. Optimal paths 
in networks have also been studied with exponential 
disorder, where the optimal path is the path between 
two sites that minimizes the total weight J2 path e aXij 
El IH 111 13 E3, where the sum is over the bonds 
(ij) along the path. The length of the optimal path £ opt 
has been shown to scale with the system size as L d °v t for 



the strong disorder limit (a — > 00) [21|. where a single 
bond dominates the optimal path (and conductance as 
we see below). The strong disorder limit only has been 
related to critical percolation 0, H H 0] . 

Here we show that for exponential disorder, the flow 
paths for all values of a are controlled by critical percola- 
tion and by the scaling properties of the optimal path in 
the strong disorder limit. Indeed, the resistance of each 
path is equal to the sum of its resistances. When a — > 00 
the resistance of each path is dominated by the largest 
resistance on this path exp(aa; m ax)- Almost all currents 
must go along the path which minimizes x max . We denote 
this min-max value of disorder as x\ = min a n paths ^max- 
Among all the paths which go through the bond with x\ , 
the maximum-current goes along the path which min- 
imizes the second largest value of disorder x' max , and 
so on. Thus the algorithm of selecting the path with 
the maximum-current is equivalent to selecting the opti- 
mal path in the strong disorder limit (ultrametric algo- 
rithm [l6j). As a — > 00 the maximum-current path coin- 
cides with the optimal path in the strong disorder limit. 
On the other hand, since all values x^ on the maximum- 
current path are below x\, this path must belong to the 
percolation backbone with concentration p equal to the 
fraction of bonds whose x^ < x\ The value of p at 
which percolation between two edges of the system does 
occur has a narrow distribution with a mean of p = y r 
and a standard deviation that scales as ~ L~ x l v |23j . 
where p c is the critical percolation threshold, L is the 
linear system size, and v is the connectedness length ex- 
ponent. Thus the value of x\ also must have a narrow 
distribution of width ~ L~ x l v . 

Next we estimate the value a at which the maximum- 
current path starts to bifurcate. Consider the paths 
which do not pass through bond x\ as if this bond 
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FIG. 2: (a) Plot of Pit \ L, u) versus £ for square lattices 
with fixed u = L/a v — 10 and different values of L. (b) Plot 
of P(£ | L, u)L ' versus l/L 1 - 22 . Three families of curves 
are shown and each family has the same ratio u = L/a v : 
u = 10 [L = 250(o), L = 300(D), L = 350(o)]; u = 1.26 
[L = 30(A), L = 40(<), L = 60(v)]; « = 0.25 [L = 15(+), 
£ = 20(x), L = 27(*)]. The distribution curves with the 
same ratio of u collapse both in weak disorder such as u — 10 
and in strong disorder such as u — 0.25, as well as in the 
intermediate regime u — 1.26 for a > 10 and L > 15. We 
compute all the data with 1000 realizations of disorder and 
10 5 tracers for each realization. 
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FIG. 1: Current maps for the same configuration of disorder 
on a 15x15 square lattice, with different values of a: (a) a = 
5, (b) a = 20, and (c) a = 45. Each bond represents one 
resistor. The dot density of each bond increases when the 
bond current increases. The source has coordinate (0,7) and 
the sink (14,7). (d) The corresponding optimal path for the 
same configuration of disorder and for a = 45. The similarity 
between (c) and (d) suggests a relation between current flow 
paths for large a and the optimal path. 



has been cut [2|. The maximum-current will then pass 
through bond X2 > x\, which is characterized by the 
same narrow distribution. Hence (x2 — xi) is of the order 
of L~ x l v . These paths become competitive with the true 
optimal path if its resistance exp(ax2) becomes of the 
same order as exp(axi) or if a(xi — x\) ~ aL~ x l v w 1. 
This condition determines the crossover from weak to 
strong disorder. If L -C a v ', the disorder is strong and 
the maximum-current path does not bifurcate. If L 3> a" 
the disorder is weak and the maximum-current path can 
bifurcate. Moreover, the value £ ~ a v determines the 
connectedness length below which the disorder is strong 
and the maximum-current path is determined by the 
unique optimal path and above which the maximum- 
current path bifurcates. 

To confirm these analytical predictions, we study the 
problem numerically. Define the electric potential at 
node i of the lattice as Vi , and set the potentials at source 
and sink as Va = 1 and Vb = 0, we numerically solve 
the set of Kirchhoff equations for all Vi |2J|. We begin 
by building an intuitive understanding of the effect of 
changing the strength of disorder on current flow. Fig- 
ures [D^a), (b), and (c) show, for different values of a, the 
magnitudes of the bond currents represented by the den- 
sity of dots on each bond. We see that the set of bonds 



carrying most of the current decreases as a increases, so 
that only a few current paths dominate. This confirms 
earlier findings that for large a, one or very few paths 
dominate the current flow 0,0 EJ I n Fig- Did), we 
plot the optimal path for the same disorder realization. 
The similarity between the path of the current carrying 
bonds in Fig. ^c) and the optimal path in Fig. ^d) ex- 
hibits how these two quantities are related in the strong 
disorder limit and supports the argument above that the 
maximum-current path coincides with the optimal path. 

Figure ^ illustrates that the paths used by the cur- 
rent are intimately related to the disorder of the system. 
Therefore, we study the ensemble of current paths on 
the lattice by performing tracer dynamics with the par- 
ticle launching algorithm |25| . For a given realization, 
all bond currents are determined by Kirchhoff equations 
and then tracers are injected into node A and extracted 
at node B. At a given node, the tracer follows the bond 
from node i to j with probability 

^ = (2) 

where j runs over all the neighbor bonds of node i, Jij = 
Iij if lij > 0, and = if Uj < 0, so that only "out" 
currents are taken into account. 

To understand the behavior of the current flow in the 
presence of disorder in all ranges of disorder, we calcu- 
late the length distribution of all tracer paths, P(l \ L,a), 
from A to B for a system of linear size L and disorder 
strength a. We first fix u = Lja v and calculate the dis- 
tribution P u (£) = P(£ | L,a) for different system sizes 
L and the corresponding values of a = {L/u) 1 / v . Wc 
obtain weak disorder when u 3> 1 and strong disorder 
when u <ti 1, as found for the optimal path in networks 
|2rJ | and as shown below for current flow. Moreover, we 
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FIG. 3: Schematic illustration of the flow path inside and 
outside the critical regimes in the weak disorder case, £ = 
a" <C L. The parameter u — L/a v determines the number of 
such units (of size f) in a linear size L. While the total length 
of the flow path is linear with L, for distance y < £ (inside 
the critical regime) , we expect £ ~ y d °P t . 



find that u is the only parameter that characterizes the 
disorder and thus determines P u (l). 

In Fig. Ota) we show three normalized distributions 
P u (£) with u = 10 (weak disorder), which collapses to 
a single curve as shown in Fig. Hfb) • Figure |3b) also 
shows two other peaked curves with u = 1.26 (close to 
the crossover) and u = 0.25 (strong disorder). Each curve 
shows the collapse of three distributions with different 
system sizes L but the same value of u. This collapse 
implies that P(£ \ L, a) is controlled by a single parame- 
ter u 



P(£\L,a) 



1 



L d °pt 



In 



(3) 



We confirmed this scaling numerically for values of u be- 
tween u = 10 (weak disorder) and u = 0.25 (strong dis- 
order) for a > 10 and L > 15 |27j. 

To understand why v and d opt play an important role 
in determining the length of the current flow path in 
weak disorder as well as in strong disorder, we suggest 
the following theoretical argument. In the weak disor- 
der regime, there is a characteristic length £ ~ a v below 
which strong disorder exists and critical percolation plays 
a crucial role [2^ . We thus expect that for length scales 
up to £ the tracers travel on strong disorder path seg- 
ments with a typical length of ~ £ d °pt, and a tracer 
length deviation of <j£ ~ £ do Pt (illustrated in Fig.^J). For 
a system of linear size L in weak disorder, the ratio of 
the system size to the connectedness length u ~ L/£ 
roughly indicates the number of independent strong dis- 
order tracer path segments within a complete tracer path 
from source to sink. The total length is obtained by mul- 
tiplying u by the length of a segment, £ rf opt. pj e g mn g £* 
as the maximum of P{£ \ L,a), we thus predict that in 



(a) 







(b) 



0.1 



1.0 



10.0 



0.1 



1.0 



10.0 



L/a 



L/a 



FIG. 4: (a) Log- log plot of t/L 1 ' 22 versus L/a" for different 
values of a: 10(o), 15(D), 20(o), 25(A), 30(<), 34(v) and L 
changes from 20 to 200. The slope of the dashed line is -0.21, 
in agreement with Eqs. © and ©. (b) The same as (a) but 
for a. The slope of the dashed line is -0.72, in agreement with 
Eqs. © and ©. 



the weak disorder I* is 

^«opt = £°opt<i/-<*opt , 



£* £ r -.£<2nnt T dr.T.t„.l — d n 



(4) 



where £ is the mean average path length of the tracers [23 . 
l3(ij |. Thus for all values of u, £* can be written in a unified 
form 



t „ L d °rtg e (u), 



(5) 



where gi{u) is a scaling function that satisfies (from 
Eq. Q) 



y} d OP t U > 1 

1 u < 1. 



(6) 



The arguments leading to Eq. Q for weak disorder, also 
imply that the standard deviation a of £ scales as 



a ~ V< dopt = L^ptu 1 / 2 -^*, 
and for all values of u 

with 



9cr(u) 



u V2-rf op t u > 1 
1 u < 1. 



(7) 
(8) 
(9) 



To test Eq. 10) we plot ^*/i d °P t as a function of u 
in Fig. Ufa). We find that the best scaling is obtained 
for d op t = 1-22, the predicted value. When L ^> a" 
(u 3> 1), 5^(u) is asymptotically a power law function 
with an exponent —0.21 ± 0.02, which is within the error 
of the predicted value 1 — <i pt = —0.22 ± 0.01 (from 
Eq. ©). Similarly, in Fig. EJb), we plot ajL d °^ as a 
function of u = L/a u and find that g a is asymptotically 
a power law with an exponent —0.72 ± 0.02 as predicted 
in Eq. 0. All these results strongly support our picture 
of critical percolation regimes of size £ ~ a v . 

Equations (J5J and J§J state that tracer path length 
scales with system size L in the same way as the optimal 



4 



path length for all values of u. For u <C 1, t ~ L d °v t 
and the path is a fractal with the same exponent <i op t 
as for the optimal path length £ opt . In weak disorder 
(u ^S> 3LL, we obtain £* ~ L as we do for self-affine struc- 
tures This is consistent with the interesting possi- 
bility that they belong to the same universality class. As 
u <C 1, current flows only along the optimal path, which 
explains the existence of the bottleneck at the percolation 
threshold p c 0,0,1111. 

Our results also explain the simulation results of 
Ref. 0] for the scaled plot log(i? cut /i?) as a function 
of a/L 1 / 13 , where R is the equivalent resistance of the 
2D random resistor lattice and i? cu t is the equivalent re- 
sistance of the system after cutting the bond with the 
maximal local current. Before cutting this bond, the 
equivalent resistance in the strong disorder limit is dom- 
inated by the maximal resistance along the optimal path 
R ~ e aPc [1 HQ. After cutting this bond, the cur- 
rent will reorganize to follow a new optimal path on 
which the dominant resistance is R ~ e ap , making the 
ratio Rcut/R ~ e a{p - p ^ = e aSp . Using the relation 
S p ~ £-V" ~ L- 1 /" [H, we find that 



R 



cut pa/L 1 '- 



R 



(10) 



This result also analytically supports our assumption 
that the ratio L/a v characterizes the disorder and de- 
termines the properties of current flow. 

In summary, we find that the tracer path length £ 
in flow in the presence of exponential disorder behaves 
similarly to the optimal path length £ op t, and even has 
the same scaling exponents (c£ opt for u <C 1 and one for 
u> 1). Moreover, we also find that when the disorder is 
weak and I ~ L, there is a connectedness length £ ~ a u , 
where strong disorder and critical percolation exist for 
regimes smaller than £. As a result, the probability dis- 
tribution of I is determined by the ratio Lja v ', which is 
the number of units of size ( ~ o" in a linear size L. 
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